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The equation of state (EOS) of dense matter and neutrino mean free path (NMFP) in a neutron 
star have been studied by using relativistic mean field models motivated by effective field theory 
(ERMF). It is found that the models predict too large proton fractions, although one of the models 
(G2) predicts an acceptable EOS. This is caused by the isovector terms. Except G2, the other two 
models predict anomalous NMFP. In order to minimize the anomaly, besides an acceptable EOS, a 
large M* is favorable. A model with large M* retains the regularity in the NMFP even for a small 
neutron fraction. 



PACS numbers: 13.15.- 



25.30.Pt, 97.60.Jd 



The finite-range (FR) (see Refs. J1|1|1I3) and point- 
coupling (PC)(see Refs. 0,IE0ilEI3) types of relativistic 
mean field (RMF) models have been quite successful to 
describe the bulk as well as single particle properties in 
a wide mass spectrum of nuclei. 

The early version of RMF-FR is based on a Lagrangian 
density which uses nucleon, sigma, omega and rho mesons 
as the degrees of freedom with additional cubic and quar- 
tic nonlincaritics of sigma meson. For example NLZ, 
NLl, NL3 and NL-SH parameter sets belong to this ver- 
sion. Recently, inspired by the effective field and den- 
sity functional theories for hadrons, a new version of this 
model (ERMF-FR) has been constructed It has 

the same terms like the previous RMF-FR but with addi- 
tional isoscalar and isovector tensor terms and nonlinear 
terms in the form of sigma, omega and rho mesons com- 
bination. One of parameter sets of this version is G2. 
Besides yielding accurate predictions in finite nuclei and 
normal nuclear matter 0,0,01, G2 has the demanding 
features like a positive value of quartic sigma meson cou- 
pling constant that leads to the existence of lower bound 
in energy spectrum of this model and to the miss- 

ing zero sound mode in the high density symmetric nu- 
clear matter Moreover, the agreement of the nuclear 
matter and the neutron matter equation of states (EOS) 
in high density of G2 with the Dirac Brueckner Hartree 
Fock (DBHF) calculation 0, [III is better than those of 
NLl, NL3 and TMl models (the standard RMF-FR plus 
a quartic omega meson interaction). 

The difference between RMF-PC and RMF-FR is due 
to the replacement of mesonic interactions in the FR 
model by density dependent interactions. It is evident 
that RMF-PC and RMF-FR serve similar quality inpre- 
dicting finite nuclei and normal nuclear matter [a, S 13 ■ 
This is due to the fact that "finite-range" effects in RMF- 
PC model are effectively absorbed by the coupling con- 
stants. Therefore in connection with different treatments 
of "finite-range" in both models, studying the behavior 
of the PC model in high density should be interesting. In 
this report, we choose the VA4 parameter set of Ref. 
(ERMF-PC) because it can be properly extrapolated to 
the high density and it has also density dependent sclf- 



and cross-interactions in the nonlinear terms. 

So far the EOS of a neutron star has not been known 
for sure However, recently ^3 the flow of mat- 

ter in heavy ion collisions has been used to determine 
the pressure of nuclear matter with a density from 2 
until 5 times the nuclear saturation density (po)- Ref- 
erence has found that these data can be explained 
only by the variational calculation of Akmal et al. |l6l |. 
Unfortunately, this interaction cannot be successfully ap- 
plied to the case of finite nuclei 0| . Reference ^Jl found 
that the EOS predicted by G2 is in a gree ment with data. 
This result is remarkable, since Ref. [13 states that the 
minimal requirement for an accurate neutrino mean free 
path (NMFP) is a correct prediction in the low density 
limit, as well as the consistency with the correspond- 
ing EOS. On the other hand, one should remember that 
many-body corrections arc important but they depend 
on the model and the a ppr oximation of strong interac- 
tion used El 111 111 iiiM El mm mm. 

According to Refs. IH IH all RMF-FR models yield 
lower threshold densities for direct URCA process than 
those of variational calculations |Tal . In the neutron star 
cooling model, Migdal et al. [23| treated this fact as a 
fragile point of RMF-FR models. So, they disregarded 
direct URCA from their analysis but Lattimer et al. 
used this fact to develop their direct URCA scenario. 

Therefore, in this report we will compare the neutron 
matter prediction in high density from the G2, NLZ and 
VA4 models in order to check the result of Ref. 0] and 
the possibly different predictions from ERMF-PC and 
ERMF-FR due to the different treatment of the "finite- 
range effects" . Furthermore, the agreement between the 
G2 EOS with experimental data has motivated us to cal- 
culate NMFP using this model for direct URCA pro- 
cess. A similar assumption as in Ref. is used, i.e., 
the ground state of the neutron star is reached once the 
temperature has fallen below a few MeV. This state is 
gradually reached from the later stages of the cooling 
phase. The system is then quite dense and cool so that 
zero temperature is valid. In this case the direct URCA 
neutrino-neutron scattering is kinematically possible for 
low energy neutrinos at and above the threshold den- 
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sity when the proton fraction exceeds 1 /9 [23 or slightly 
larger if muons are present. Furthermore, the absorp- 
tion reaction is suppressed. For simplicity, we neglect 
the RPA correlations. 

The effects of self- and cross- interactions terms and the 
treatment of finite-range in high density can be observed 
by extrapolating the EOS which is presented by the neu- 
tron matter pressure P and the effective mass M*, as 
shown in Figs. ^ a-nd El where we compare the results 
obtained from the G2 g^, NLZ and VA4 models 
as a function of ps/ Po- 



ll is found that the nuclear matter EOS of VA4 is stiffer 
than those of NLZ and G2, even for pB less than po- How- 
ever, the G2 EOS is softer than the NLZ one at the high 
density but not at the low density. This fact emphasizes 
the result of Ref. jlj that the crucial role of self- and 
cross-interactions of meson exchange model is to soften 
the EOS at the high density. 

It is shown in Fig.|2]that for 1 < ps/po < 5, the effec- 
tive mass > Myj^^, but ioi pB / po > 5 one observes 
that < ^vAi- This indicates that quantitatively 

M* depends on the model. We note here that the effec- 
tive masses of G2 and VA4 depend on self- and cross- 
interaction terms implicitly. We also note that other 
mechanisms could also produce a larger M* , e.g., in the 
Zimanyi-Moszkowski and linear Hartree-Fock Walecka 
models |22j . where those terms are not present. Although 
those models give a regular NMFP, they are quite unsuc- 
cessful in finite nuclei applications, especially in predict- 
ing the single particle spectra of nuclei 0|. Therefore, 
it is interesting to check whether or not the relation be- 
tween a large M* and a regular NMFP also appears in 
the case of ERMF models. 

Now, we calculate the NMFP of the neutron star mat- 
ter by ernploying G2, VA4 and NLZ models. Following 
Refs. I21II23I, we start with the neutrino differential scat- 
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FIG. 1; Equation of states (EOS) of the neutron matter. 




FIG. 2: Effective masses (M*) of the neutron matter. 
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Here E^^ and E"^ are the initial and final neutrino energies, 
respectively, Gp = 1.023x10~^/M2 is the weak coupling, 
and M is the nucleon mass. The neutrino tensor L^^, can 
be written as 



(2) 



where k is the initial neutrino four-momentum and q ~ 
{qo,q) is the four-momentum transfer. The polarization 
tensor 11^'', which defines the target particle species, can 
be written as 
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where j= n,p,e ,/i . G{p) is the target particle prop- 
agator and p = (j)o,P) is the corresponding initial four- 
momentum. The currents are "i^{Gy — C^75). The 

explicit forms of G^ (p), Cy and C;^ of every constituent 
and also their explanations can be found in Ref. [2^ . 
The NMFP (symbolized by A) as a function of the ini- 
tial neutrino energy at a certain density is obtained by 
integrating the cross section over the time- and vector- 
component of the neutrino momentum transfer. As a 
result we obtain [2lll23 
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Since in our study we assume that the neutron star 
matter consists only of neutrons, protons, electrons, and 
muons. the relative fraction of each constituent should 
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FIG. 3: Neutron fraction in the neutron star matter with and 
without isovector terms. 



FIG. 4: Neutrino mean free path (NMFP) in the neutron star 
matter. 



be taken into account in the NMFP calculation. The 
relative fraction is determined by the chemical potential 
equilibrium and the charge neutrality of the neutron star 
at zero temperature. The neutron fractions for all models 
are shown in Fig. |31 

Qualitatively, all parameter sets have similar trend in 
fraction of each constituent, i.e., when the neutron frac- 
tion is decreasing, other constituent (p, e~,/x~) fractions 
are increasing. Quantitatively, isovector terms are re- 
sponsible for the high proton fraction. G2 has smaller 
neutron fraction than VA4 and NLZ. Therefore, even 
though G2 has an acceptable EOS, it has a too large 
proton fraction. This fact leads to such a low threshold 
density for direct URCA process. We note that this fact 
is ruled out by the analysis of the neutron stars cool- 
ing data [26l Isil Is^ . Thus, this result indicates that 
significant improvements in the treatment of isovector 
sector of ERMF-FR are urgently required. Variational 
calculation of Akmal et al. jig allows for a direct URCA 
process only for pb/po > 5. Linear Walecka (hnear FR) 
and Zimanyi-Moszkowski (derivative coupling) Hartree- 
Fock models of Ref. yield a higher critical density 
for the direct URCA process. Isovector contributions of 
these models do not drastically change the proton frac- 
tion. But on the other hand, all Hartree-Fock models 
of Ref. |23| are unable to give a good prediction in finite 
nuclei, especially in the single particle properties 33., 34] . 
It may be interesting to see also the consistency of their 
EOS with experimental data [l5l |. 

The NMFP for all models can be seen in Fig. ^ Here 
we use a neutrino energy of i?^ = 5 MeV. In general, 
from medium to high density, Xnlz is larger than Ag2 
and XvAi- In high density region we clearly see that 
Ag2 ~ XvAi- The NMFP difference among all models 
appears to be significant around 1 < pb/ Po ^ ^ (medium 
density). For pb/po smaller than 1, Xva4 ~ Ag2 Xnlz- 
In other words, in the limit of low density all parameter 
sets serve similar A prediction as we expected. 
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FIG. 5: Neutrino mean free paths (NMFP) in the neutron 
star matter as functions of M* and neutron fraction for all 
parameter sets. 



In Fig. [3 we show the dependence of A with respect 
to the M* and neutron fraction. Obviously, NLZ has a 
maximum NMFP at M* w 200 MeV and neutron frac- 
tion w 0.75. These lead to a bump in Xnlz as shown 
in Fig.^J On the other hand, G2 demonstrates no maxi- 
mum in M* and neutron fraction dependences, leading to 
a smoothly decreasing function of Ag2 displayed in Fig.^ 
For comparison, previous NMFP calculations by using 
all Hartree-Fock models 0| showed also no anomaly. In 
these models, the predicted NMFP falls off faster than 
that of the Hartree type model as the density increases. 

In conclusion, the EOS and NMFP of ERMF mod- 
els in the high density states have been studied. It is 
found that the ERMF-FR and ERMF-PC models have 
different behaviors in high density and even by using a 
parameter set that predicts an acceptable EOS, the cal- 
culated proton fraction in neutron star is still too large. 
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Isovector terms are responsible for this. Therefore, im- 
provements in the treatment of the isovector sector of 
ERMF-FR should be done. Different from the Hartree- 
Fock calculation of Ref . , only the parameter set with 
an acceptable EOS (G2) has a regular NMFP. In order to 
minimize the anomalous behavior of A, a relatively large 
M* in RMF models is more favorable. It seems that the 



relatively large M* in the ERMF models at high den- 
sity originates from the presence of the self- and cross- 
interactions in nonlinear terms. The RMF models with 
relatively large Af * retain their regularities partly or fully 
even for a small neutron fraction. 
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